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Abstract 

We report on the time dependent solutions of the g—generalized Schrodinger equation pro¬ 
posed by Nobre et al. [Rhys. Rev. Lett. 106, 140601 (2011)]. Here we investigate the case of 
two free particles and also the case where two particles were subjected to a Moshinsky-like 
potential with time dependent coefficients. We work out analytical and numerical solutions 
for different values of the parameter q and also show that the usual Schrodinger equation 
is recovered in the limit of g ^ 1. An intriguing behavior was observed for g = 2, where 
the wave function displays a ring-like shape, indicating a bind behavior of the particles. 
Differently from the results previously reported for the case of one particle, frozen states 
appear only for special combinations of the wave function parameters in case of g = 3. 


Highlights 

Solutions of a nonlinear Schrbdinger equation; 
Time dependent wave functions; 

Free particles and Moshinsky-like potential. 

Email address: lgaalvesOdfi.uem.br (L. G. A. Alves) 


Preprint submitted to Physica A. 


February 26, 2015 



1. Introduction 


Tsallis statistics nisi. based on the entropy 



( 1 ) 


where fc is a constant and g is a real parameter, has shown to be useful in the description of 


several systems where the additivity fails such as systems characterized by non-Markovian 
processes lailiiiiT!, nonergodic dynamics and long-range many-body interactions [5] . The 
nonadditive characteristic of Sg EEiiini has also brought new insights to several works 
in the area of complex systems pun] and contributed to diverse problems of quantum 
mechanics [iiiiiaiisiiiiiiiaiisiiniiiiiiisiEQiEi]. Recently, Nobre, Rego-Monteiro and 
Tsallis have proposed a nonlinear Schrodinger equation [22| (NRT equation) that admits 
soliton-like solutions (g-plane waves) with possible applications in several areas of physics, 
including nonlinear optics, superconductivity, plasma physics, and darkmatter. Since the 
applicability of linear equations in physics is usually restricted to idealized systems [22| , the 
advance in the understanding of real systems makes necessary the use of nonlinear equations 
several times [23] • For a system with particle of mass m, the NRT equation can be written 


as 



1 ^2 



( 2 ) 


2 — q 2m 


where the scaling constant $o guarantees the appropriate units for the different physical 
terms appearing in the equation, i is the imaginary unity and h is the Planck’s constant. 
An interesting aspect related to this equation is the non-Markovian effect emerging from 
the presence of the nonlinearity, i.e., q ^ 1. A similar situation is found in the porous 
media equation [23], whose solutions may exhibit a short or a long tailed behavior and 
are connected to anomalous diffusion [23]. This intriguing equation has been the focus of 
intensive research in recent years [261I2ZIE811291I30]. For instance, a family of time dependent 
wave packet solutions have been also investigated in Ref. ra. a quasi-stationary solution for 
the Moshinsky model has been proposed in Ref. [2H] , and a classical field-theoretic approach 
has been considered in Refs. [29] [30] . 


2 







Here we study the problem of two particles in the context of the NRT equation focusing on 
time dependent solutions. We worked on the case for free particles (Section 2 ) and particles 
subjected to a time dependent Moshinsky-like potential (Section 3). For these situations, 
we observe that the nonlinearity creates an entanglement between the particles, which is not 
present in the usual scenario and it is essential for describing physical reality [HI] inherent 
to quantum mechanics [32l |33| . We summarize our results and conclusions in Section 4. 


2. Free particles 


Let us start our discussion by considering the g-generalized Schrbdinger equation for 
a system of two particles in the absence of potential by assuming mi = m 2 = m and 
defining = ^/'Fo without loss of generality. The time dependent Moshinsky-like potential 
is discussed in the next section. For this case the NRT equation takes the following form 


d 






92 92 


'lp{Xi,X2,t) 


2-q 


(3) 


where Xi and X 2 represent the particles coordinates. 

Similarly to the situation worked out in Ref. [22], Eq.([^ has the g-plane wave as solution 


'ilj{xi,X2,t) = expg [i{kixi - (uR)] exp^ [i{k2X2 - U2t )], (4) 

where is the g-product [HI]. This result shows that the nonlinearity of the Eq.Q en¬ 
tangled the particles, leading us to different behaviors dependent on q (see Eig.([^). It is 
worth noting that obtained from the NRT equation does not has a probabilistic inter¬ 
pretation for g 7 ^ 1. In order to recover that probabilistic interpretation, we could follow 
the procedures proposed by Rego-Monteiro and Nobre [22], and introduce a second field 
defined by means of an additional non-linear equation coupled with the NRT equation for 
g 7 ^ 1. Eurthermore, recovering the probabilistic interpretation of the NRT equation under 
the action of an spatial and time dependent potential (such as the Moshinsky-like poten¬ 
tial worked out in the next section), in contrast with the free particle case, still is an open 
problem. Even in the free particle case, to obtain general time-dependent solutions for the 

NRT equation represents a cumbersome task due to the intrinsic complexity of xp that still 
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Figure 1: Behavior of Eq. @ for two different values of q in order to illustrate the entanglement of the 
particles due to the nonlinearity of the NRT equation (Eq.Q). 


need to be incorporated in this second equation for obtaining the second field. For this 
reason, we have focused our study only on the solutions for the NRT equation, that is, on 
'll;. Although this simplified approach lacks an exact probabilistic interpretation, for q not 
so far from one we could expected that will keep some information on the probability 
of find particles in a given position and time; also, this analysis can provides new insights 
and a better understating of the NRT equation. 

In order to explicit solve Eq.Q, we use the g-Gaussian package as an ansatz, 

'ip{xi,X 2 , t) = [l - (1 - q) [a{t) x\ + h{t) xl + c{t) xi X 2 + d(t))] (5) 


where a, b, c, and d are appropriate time dependent coefficients to be found. The presence of 
a crossing term in Eq. (|^ covers a more general scenario characterized by an initial situation 
of mixing between the variables xi and X 2 . 

By substituting ([^ in (|^, the left and right sides of the NRT equation becomes, 
d 

ih—ip{xi,X 2 ,t) = —ih (x^a'(t) + x\b'{t) + Xi X 2 c\t) + d\t)) 'ij;{xi,X 2 ,ty, (6) 


and 


1 


d‘^ d‘^ 

+ 


Y 




'lp{xi,X2,tY —V'(xi,X2,t)'^ ) = 


2 — q 2m \ dx\ dx\ 

2a{t) ((g - l)b{t) (xf + x^) + (g - 3)a:iX2c(t) + (g - l)d(t) + l) + 2(g - 3)x‘la{tY+ 


2b(t){{q — 3)xiX2c{t) + (g — l)d{t) + 1) + 2(g — 3)xlb{tY — c(t)^ + xl). (7) 
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Thus, Eq.([^ satisfies the NRT equation with a, b, c and d being solution of the following 
set of coupled ordinary differential equations: 

ih 


a'{t) = -— (2 a(t)((g - 3) a{t) + {q - 1) b(t)) - c{tf) , 


2 m 
ifi 


b’(t) = — (2 bmg-l)a(t)+ (9-3) !>(())-c(tf). 


2 m 
jh 

c'it) = _(g-3)c(t)(a(t) + 6(t)), 
m 
jh 

d'{t) = -{a{t)+b{t)){{q-l)d{t) + l). 

m 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Note that there is a relationship between these parameters. By considering the difference 
a'{t) — b'{t), we obtain 


a'it)-b'it) = 


c'{t) 

c{t) 


{a{t)-b{t)). 


Thus, we can write the following relationship for a, 6 and c: 


a{t) = b{t) + K c(t), 


( 12 ) 


(13) 


where k = (a(0) — 6(0 ))/c(0). 

A general analytical solution for this set of equations with arbitrary initial conditions is 
a hard task; however, note that any initial state with c(0) = 0 lead us to c{t) =constant 
for all later times. Likewise, any initial state with o(0) = 6(0) will preserve this symmetry 
for all later times (see Eq.([I^). Also, if c(0) = 0 and a(0) = 6(0), we have the following 
solutions: 

'ihcQt 


Co 


a(t) = - , _ 

b{t) = a{t), 
c{t) = Co, 


tanh 


m 


■ \/Q ~ + 2 Co Cl q — 2 j , 


d{t) = 


1 


1-q 


C 2 e 


2ih (g-1) t 


COS 


hcot 


m 


~ 2 + 2 i Co Cl q — 2 


g-1 

2(g-2) 


(14) 

(15) 

(16) 

(17) 


where, co, ci and C 2 are constant that depend on the initial conditions. Naturally, it is 
possible to obtain other numerical and analytical solutions for a, 6, c and d to investigate 
the behavior of \'tp{xi,X 2 ,t)\'^ for the cases that we shall present. 


















2.1. Limit case q ^ 1 


Let us discuss the case that recovers the linear Schrodinger equation for a systems of two 
entangled free particles. In this case, the set of differential eqnations for o, b, c and d is 

a'it) = + £©!), ( 18 ) 

m = ( 19 ) 

2 if) 

c'{t) = - c{t) {a{t) +b{t)), ( 20 ) 

Tin 

ifi 

d\t) = - {a{t) + b{t)). (21) 

im 

Considering the difference a'{t) — b'{t), Eq.( [^, and using this result in the above set of 
equations, we can rewrite c'{t) as 

ih 

c'{t) = - (4 b{t) c{t) + 2 K c{t)^) . ( 22 ) 

If we sum b'{t) and c'{t) we have that, 

b'{t) + c'{t) = —^ + 2 6 (t) c{t) + ■ (23) 

Notice that k is a free parameter that depends only on the initial conditions of a, b and c. 
For simplicity, we choose k = 3/4, yielding 

O jf) 

b\t) + c\t) = -^( 6 (t) + c(t))' (24) 

m 

and, therefore 


b{t) 


m 

2iht 


c{t). 


We solve the set of eqnations for the parameters a, b, c and d, leading to: 


a{t) 

b{t) 

c{t) 

d{t) 


2mh + i \ m? t 

—2 Xmt^h + lliU^ T 
m 2m 

2ifit t{2Xmt — llzh)’ 
2 m 

t{2Xmt — llih)' 
llh, + AiXmt 

22th + 4:i Xmt"^’ 
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(26) 

(26) 

(27) 

(28) 
(29) 












where A depends on the initial conditions. It is worth noting that when t —?• cx), we have 
a{t) ~ b{t) ~ d{t) ~ 1/t and c{t) ~ 1/t^. We have also evaluated a numerical solution by 
considering the initial conditions o(0) = 1, 6(0) = 1, c(0) = 1, d{0) = 1 and plotted the 
evolution of | xp(xi,X 2 ,t) p for different times in Fig.Q. In this case, the contribution of 
the coupling term, c{t), vanishes faster than the others for long times, making the shape of 
the wave function similar to the case of one free particle rni. 



Figure 2: Evolution of the | '0(xi,X2,t) p for the entangled system of free particles in the linear limit case, 
q = 1. Here, we have used m = 1, = 1, and the following initial conditions: a(0) = 1, 6(0) = 1, c(0) = 1 
and d{0) = 1. We have employed the In | '0(xi,X2,t) p in order to smooth the surfaces. Figure (A) shows a 
3D-plot and figure (B) shows a density-plot of the In | '0(xi,X2,t) p in function of the positions xi and X 2 
of the two particles. 

2.2. Evolution of the wave funetion for q=2 

For the nonlinear case g = 2, the set of differential equations takes the form: 


a\t) = 

1 h 

Zm 

(30) 

h'{t) = 

1 n 

9^(2K^)(«6) Kt)) 

Zm 

(31) 

c'it) = 

ih 

- c{t) {a{t) + b{t)), 

m 

(32) 

d'{t) = 

if) 

-{d{t) + l) {a{t) + b{t)). 
m 

(33) 
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Because the inherent complexity of this set of equations, we assume that a{t) = b{t) in order 
to simplify the problem. Under this assumption, we have the following solutions: 


a{t) 

c{t) 

d{t) 


b{t) = ± J^tan ( ^/^C 2 - 


m 


=F\/ 2 ci sec ( \/2 ci C 2 — 

\/2 Cl 


i\/2 Cith 


-1 + C 3 cosh 


m 


m 

{th + ic2m) ) , 


(34) 

(35) 

(36) 


where ci, C 2 and C 3 are constants that depend on the initial conditions. In particular, 
the constant Ci can be determined via Ci = 2 a^( 0 ) — c^( 0)/2 and the others by using the 
previous results. In Fig.([^, we show a numerical solution for | tp{x, y, t) p where we observe 
the emergence of two peaks around t = 50, which evolve to a ring-like shape mt = 100 (see 
Fig.([^)) and bind the xp function of the two particles. This behavior was also found in the 
case of one particle [ 12 ]. 



X2 


X2 


X2 


Figure 3: Evolution of the | ip{xi,X 2 ,t) p for the entangled system of free particles in the case q = 2. Here, 
we have used m = 1, h = 1, and the following initial conditions: a(0) = 1, 6(0) = 1, c(0) = 1 and d(0) = 1. 
We have employed the ln| ■ip{xi,X 2 ,t) p in order to smooth the surfaces. Figure (A) shows a 3D-plot and 
figure (B) shows a density-plot of the In | 'ip{xi,X 2 ,t) p in function of the positions Xi and X 2 of the two 
particles. 
























2.3. Analytical solution for q=3 and the absence of “frozen” solutions 

For the case of g = 3, the set of differential equations for the parameters o, b, c and d 
can be reduced to 


aft) 

= 2 m 

(37) 

bft) 

Zm ^ 

(38) 

eft) 

= 0 , 

(39) 

dft) 

= —( 2 d(t) + l)(a(t) + 6 (t)). 

TTh 

(40) 


Notice that aft) = b'{t) and the solution for c'ft) is straightforward. Also, we have that 
aft) = bft) by imposing the same initial conditions for Eqs.(33) and (34). Under this 
assumption, the time dependent parameters are given by 


aft) 

c{t) 

dft) 


b{t) 


Co fl - 

2 VI + J ’ 


Co, 

-il + 2do) 

2 2 


(41) 

(42) 

(43) 


where oq, cq and do depend on the initial conditions, A = {cq/2 + ao)/(co /2 — oq) and 
a = 2h/m (see Fig. |^. Curilef et al. na have showed the existence of frozen solutions for 
g = 3 in the case of one particle. In our case, there are no frozen states for 'tp{xi,X 2 ,t) 
due to the interactions between the particles caused by the nonlinearity of the Eq.Q. We 
illustrate this aspect by plotting the evolution of the parameters o, 6 , c and d (in absolute 
values) for g = 1, g = 2 and g = 3, as shown in Eig.([^. The frozen states appears only in 
special cases where 4a(t) bft) = cft)"^ and dft) = 0. We further observe that if ao = —Co/2, 
then aft) = bft) = —co /2 and dft) ~ g-*acot^ which implies in a quasi-stationary behavior 
for '^(xi,X 2 ,t) . In this scenario, the solution “pulsates” with a period ^ and for do = 0 
we obtain another frozen state, similarly to the case of one free particle da. 


3. Moshinsky-like potential 

Entanglement features of atomic systems are often studied by considering the Moshinsky 

model [HSl ES] • This model has been intensively investigated in the recent years as a test 
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X2 


X2 


X2 


Figure 4: Evolution of the | '0(xi,X2,t) p for the entangled system of free particles in the case g = 3. Here, 
we have used m = 1, = 1, uq = 1, cq = 1 and do = 1. We have employed the In | '0(xi, X 2 , t) p in order to 
smooth the surfaces. Figure (A) shows a 3D-plot and figure (B) shows a density-plot of the In | '0(xi, X 2 , t) p 
in function of the positions xi and X 2 of the two particles. 


A q='' 



B 



c 



a(t) - b(t) - c(t) - d(t) - 


Figure 5: Evolution of the time dependent parameters. The plots (A) g = 1 and (B) g = 2 are numerical 
solutions, while (C) is the analytical solutions found for the case g = 3, Eqs. (37), (38) and (39). Note that 
a{t) = b{t) due to the choice of the initial conditions. 

for validating several approximations in atomic and molecular physics. In this context, we 
investigate the solutions of the NRT equation for a system of two particles subjected to a 
Moshinsky-like potential, i.e.. 


V{xi,X 2 ,t) = a{t) x\ + j3{t) x\ + 7(t) xi X 2 + r){t) 
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(44) 























where a, P, ^ e rj are time dependent parameters. 

Considering the potential (44) and by applying the same procedures of Section 2 , the 


ansatz ([^ is solution of the NRT equation if the time dependent parameters satisfy the 
following set of equations: 


a'{t) = 

2^^ 3)a(t)-F(g l)b{t)) c(tf)+2ma(t)), 

(46) 

h'{t) = 

C-r (ft" (2 b(t)((q - 1) a(t) + (, - 3) b(t)) - c(tf) + 2 m/3(i)) . 

Zi TTl n 

(46) 

c'(t) = 

{{q 3)h^ c{t) {a{t)+ b{t)) + m^{t)) , 

mh 

(47) 

d\t) = 

{b? {a{t) -\- b{t)) {{q 1) d{t) + 1 ) +mr]{t)) . 

Tin n 

(48) 


Naturally, we recover the free particle case ifa = l3 = ^ = r) = 0. We study numerical 
solutions for this set of equations since analytical approaches represent a cumbersome prob¬ 
lem. We have focused on numerical solutions for 'tp{xi,X 2 ,t) considering particular values 
of q, a{t) = 15(t) = r]{t) = 1 and an exponential decay for the coupling factor, this is, 
7 (t) = e~*. This latter hypothesis makes the coupling between the particles strong for small 
times and weak for long times. 

3.1. Limit case q ^ 1 

For g = 1 we recover the linear Schrodinger equation and the set of differential equations 
becomes: 


a'(t) = 

2 ^ ^ ( 4a(t)2 c{tf) + 2 m a{t)) , 

(49) 

b'{t) = 

c(t) 2 ) +2m^(t)) , 

(50) 

d{t) = 

2n^c{t){a{t) , 

(51) 

d\t) = 

^ ^ (^^(o(^) + b{t)) + m ri{t)) . 

(52) 


The solution for this case is shown in Fig. where we note that the shape of the wave 
functions are visually similar to the case of two free particles (see Fig.([^). 
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Figure 6: Evolution of the | tp{xi,X 2 ,t) ^ for the entangled system subject to a Monshinsky-like potential 
in the limit case of </ = 1. Here, we have used m = 1, h = 1, a(0) = 1, 6(0) = 1, c(0) = 1, (i(0) = 1, and 
the following functions for the Moshinsky-like potential: a{t) = 1, j3{t) = 1, •y{t) = e“*, ri{t) = 1. We have 
employed the In | ip{xi,X 2 ,t) in order to smooth the surfaces. Figure (A) shows a 3D-plot and figure (B) 
shows a density-plot of the In | ip{xi,X 2 ,t) p in function of the positions Xx and X 2 of the two particles. 

3.2. Evolution of the wave function for q=2 

In this section, we show a numerical solution for the nonlinear Schrodinger equation with 
q = 2. In this case, the set of differential equations are: 


aft) = 

— (h^ {2a{t){b{t) — a{t)) — c(t)^) -|- 2ma{f)) , 

^ f f LfL 

(53) 

hft) = 

2mh + 2m/5(t)) 

(54) 

eft) = 

^ {-h^c{t){a{t) + b{t) + nryf))) , 

(55) 

dft) = 

—— (h^i^dit) l)((2(t) -|- b(t)) mriif)') . 

(56) 


and the numerical solutions for if{xi,X 2 , t) are shown in Fig.([^. due to the interaction with 
the potential, the binding behavior appears in different times, see for example t = 50. We 
further note that the ring-like form is deformed and that | 'ip{xi,X 2 ,t) p exhibit some peaks. 
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Figure 7: Evolution of the | ' 0 (xi,X 2 ,t) p for the entangled systems subject to a Monshinsky-like potential 
for q = 2. Here, we have used m = 1, = 1, a(0) = 1, b{0) = 1, c(0) = 1, d{0) = 1, and the following 

functions for the Moshinsky-like potential: a{t) = 1, P{t) = 1, 7 (t) = 77 (t) = 1. We have employed 

the In I ' 0 (xi,X 2 ,t) p in order to smooth the surfaces. Figure (A) shows a 3D-plot and figure (B) shows a 
density-plot of the In | ' 0 (xi,X 2 , t) p in function of the positions xi and X 2 of the two particles. 


3.3. Evolution of the wave funetion for q=3 under the Moshinsky-like potential and the 
absenee of ^frozen^^ solutions 

We also investigate the existence of frozen solutions considering g = 3 under Moshinsky- 


like potential (44). For this case, the set of differential equations are: 

a'{t) = — ih? (4 a{t) h{t) — c{tY) -h 2 m a{t)) , 

2mh 

b'{t) = {4a{t)b{t) - c{tf)+2m/3{t)) , 

i 

c'{t) = 


h 


d!it) — —'z (2 6?(t) + 1) ~\~ b(t)^ + Tfi . 

mn 


(57) 

(58) 

(59) 

(60) 


Note that all the differential equations are time dependent. Even for the particular case 
where 4a{t)b{t) = c(t)^ and d{t) = 0 the parameters are time dependent, indicating that 
the Moshinsky-like potential prohibits the existence of frozen solutions. Figure [^illustrates 

\xp\^ for this case and Fig. [^displays the evolution of the parameters a, b, c and d for g = 1, 
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q = 2 and q = 3. Notice that for g = 3, the parameters approach constant values over time, 
achieving a stationary behavior. 



X2 X2 X2 


Figure 8 : Evolution of the | 'b(xi,X 2 ,t) p for the entangled systems subject to a Monshinsky-like potential 
for q = 3. Here, we have used m = 1, = 1, a(0) = 1, b{0) = 1, c(0) = 1, d{0) = 1, and the following 

functions for the Moshinsky-like potential: a{t) = 1, P{t) = 1, 7 (t) = 77 (t) = 1. We have employed 

the In I p in order to smooth the surfaces. Figure (A) shows a 3D-plot and figure (B) shows a 

density-plot of the In | 'b(xi,X 2 , t) p in function of the positions xi and X 2 of the two particles. 


A '!=■' 



Time, t 


B '1=2 



C '1=2 



a(t) - b(t) - c(t) - d(t) - 


Figure 9; Evolution of the time dependent parameters. The plots (A) q = 1 and (B) q = 2 (C) q = 3 are 
the numerical solutions. Note that a{t) = b{t) due to the choice of the initial conditions. 
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4. Summary and conclusions 


We found time dependent solutions of the nonlinear Schrodinger equation proposed by 
Nobre et al. in [Phys. Rev. Lett. 106, 140601 (2011)]. These solutions are obtained in terms 
of g-exponentials. We have investigated the case of free particles and also the case where 
the particles were subjected to a Moshinsky-like potential with time dependent coelRcients. 
We showed that the solution of the usual Schrodinger equation is recovered for g —?• 1 . We 
also investigated solutions for g 7 ^ 1 where we observe a bind behavior of for the case g = 2 
and the absence of frozen states for g = 3, differently from the results previously reported 
for the case of one particle. 
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